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ABSTRACT 

We study the escape of Lyman continuum radiation from the disks of dwarf starburst galaxies, with and 
without supershells, by solving the radiation transfer problem of stellar radiation through them. We model 
disks with Mj = 10^-10'" Mg, with exponential surface density profiles as a function of redshift, and model 
the effects of repeated supernova explosions driving supershells out of the disks, using the hydrodynamic 
simulation code ZEUS-3D. The amount of star formation is assumed proportional to mass above some 
density threshold. We vary the threshold to explore the range of star formation efficiencies,/* = 0.006, 0.06, 
and 0.6. We find that the interstellar gas swept up in dense supershells can effectively trap the ionizing 
photons before the supershells blow out of the disks. The blowouts then create galactic outflows, chimneys 
that allow the photons to escape directly to the intergalactic medium. Our results are consistent with escape 
fractions of less than 0.1 measured in local dwarf starburst galaxies, because they are likely observed while 
the starbursts are young, before blowout. We suggest that high-redshift dwarf starburst galaxies may make a 
substantial contribution to the UV background radiation with total escape fractions > 0.2, as expected if star 
formation efficiencies are >0.06. 

Subject headings: galaxies; dwarf — galaxies: formation — galaxies: starburst — H ii regions — 
ISM: bubbles 



1. INTRODUCTION 

The standard big bang theory predicts that the primordial 
gas recombines around z ~ 1 100. However, the absence of a 
Gunn-Peterson trough in quasar spectra (Gunn & Peterson 
^ 1965; Becker et al. 2001; Djorgovski et al. 2001) indicates 
CUhat the intergalactic medium (IGM) is highly ionized by 
Q redshift z ~ 5 by an ionizing background composed of 
;_i Lyman continuum radiation from quasars (accretion- 
powered massive black holes) and galaxies (young massive 
stars). The relative mixture of these two components as a 
function of redshift is poorly constrained by observations. 
. 1^ It is important to understand the nature of the background 
K> radiation and its cosmic history because most baryons 
reside in the IGM, where stars and galaxies are born. 

The contribution of ultraviolet (UV) background radia- 
tion from quasars has been studied in a number of papers 
(e.g., Madau, Haardt, & Rees 1999; Donahue & ShuU 1987; 
Shapiro & Giroux 1987; Bajtlik, Duncan, & Ostriker 1988), 
but quasars alone appear inadequate to produce enough 
UV radiation to reionize the IGM at high redshift because 
their number density declines steeply above z = 3 (e.g.. Fan 
et al. 2001). On the other hand, a large number of star- 
forming galaxies at 2 < z < 4 have recently been discovered 
using the Lyman break technique (Steidel et al. 1996), and 
the inferred star formation rate is found to stay constant 
from z = 1 up to at least ~4 (Steidel et al. 1999). This implies 
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that the contribution from galaxies increases in importance 
with redshift and must have dominated when the first stars 
and galaxies were forming (see, however, Haiman & Loeb 
1998 for the possible roles of obscured or low-luminosity 
high-redshift quasars). It is crucial to understand the 
fraction of UV photons that can escape from galaxies. 

Determining the escape fraction of ionizing photons from 
galaxies is no easy task, so there are only a limited number 
of observational estimates available. The escape fractions of 
Lyman continuum photons in four nearby starburst gal- 
axies were first measured by Leitherer et al. (1995) using the 
Hopkins Ultraviolet Telescope. They observed flux around 
900 A, which provides a reliable estimate of the total 
number of Lyman continuum photons, even if starburst 
parameters are uncertain. The sample was reevaluated by 
Hurwitz, Jelinsky, & Dixon (1997), who found upper limits 
on the escape fraction of ~0.05-0.1. The only exception is 
Mrk 66, which had an upper limit of ~0.57. The sample was 
biased toward having the most favorable conditions for a 
high escape fraction of UV photons (e.g., high Ha and UV 
luminosities but low infrared luminosities). Heckman et al. 
(2001) used the Far-Ultraviolet Spectroscopic Explorer 
(FUSE; Moos et al. 2000) to study the far-UV spectra of six 
nearby galaxies, three of which are dwarf starburst galaxies. 
They found that the strong C ii A1036 interstellar absorp- 
tion line is essentially black in all the galaxies. Since the 
opacity of the neutral intersteUar medium (ISM) below the 
Lyman edge will be significantly larger than in the C ii line, 
they estimated the escape fractions of Lyman continuum 
photons to be less than 0.06. Deharveng et al. (2001) 
recently constrained the escape fraction of Lyman contin- 
uum radiation from a starburst galaxy, Mrk 54, to be less 
than 0.062, using FUSE and foUowing the same method as 
Leitherer et al. (1995) and Hurwitz et al. (1997). 

Steidel, Pettini, & Adelberger (2001) studied Lyman con- 
tinuum emission from the composite spectrum of 29 Lyman 
break galaxies (LBGs) at z ~ 3.4. The ratio of emergent flux 
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density at rest-frame 1 500 A to that in the Lyman continuum, 
L(1500)/L(900) = 4.6 ± 1.0, implies that the escape fraction 
is very high, greater than 0.5. The ratios L(1500)/L(900) cal- 
culated for the nearby starburst galaxies still suggest that the 
escape is less easy from nearby galaxies than from the LBGs 
(Deharveng et al. 2001). The estimated comoving emissivity 
from LBGs, (1.2 ± 0.3) x 10^^ h ergs s"' Hz"! Mpc'^ 
(Einstein-de Sitter universe with the Hubble constant 
Hf) = 100 h km s"' Mpc"'), exceeds that from QSOs by a fac- 
tor of about 5 at z ~ 3. However, their sample is drawn from 
the bluest quartile of intrinsic far-UV colors, which suggests 
that it may be biased toward large escape fractions and UV 
luminosities. Very recently, Fernandez-Soto, Lanzetta, & 
Chen (2003) studied the UV escape fractions in high-redshift 
galaxies by using Hubble Space Telescope images and found 
that no more than 0.1 of Lyman continuum photons escape 
from high-redshift galaxies. 

The observations seem to suggest that only a small frac- 
tion of ionizing radiation can escape from galaxies, and 
there is a possibility that more ionizing radiation escapes 
from galaxies at higher redshift. This result appears to con- 
tradict the theoretical predictions on the escape fractions 
from starburst galaxies (Ricotti & ShuU 2000; Wood & Loeb 
2000). Ricotti & Shull (2000) studied the escape fractions of 
ionizing radiation from small spherical high-redshift gal- 
axies with dark matter masses of Mdm = 10^-10** Mr., and 
found the escape fractions to be less than 0. 1 in halos with 
AfoM^lO^ Mq. Wood & Loeb (2000) studied both low- 
and high-redshift galaxies and miniquasars with well- 
formed disks having Mdm = 10^-10^^ M^.,. Their model 
galaxies have ionizing sources distributed throughout the 
disks as ocn or rp-, and their model miniquasars have single 
ionizing sources at the center of the disks; the miniquasars 
have higher total ionizing luminosities than the galaxies. 
The escape fractions of Lyman continuum photons are 
found to be ~ 1 at low redshift for all the galaxies and mini- 
quasars and decrease to less than 0.01 for the galaxies but 
only to ~0.3 for the miniquasars at z ~ 10. Therefore, they 
suggest that the reionization must have been dominated by 
miniquasars. The study by Benson et al. (2001) supports the 
above results. 

The escape fraction of ionizing photons from our Galaxy 
studied by Dove & Shull (1994) and Dove, Shull, & Ferrara 
(2000) is consistent with that required to sustain the Reynolds 
layer. Dove & Shull (1994) studied the escape fraction from a 
collection of OB associations distributed throughout the 
Galactic disk and found it to be ~0. 1 . The studies by Ciardi, 
Bianchi, & Ferrara (2002) and Clarke & Oey (2002) confirm 
the above result. Dove et al. (2000) extended the study to 
include the dynamics of bubbles and H ii chimneys created 
by OB associations and found that the shells of expanding 
superbubbles quickly trap the ionizing flux. The estimated 
escape fraction decreases roughly by a factor of 2. 

Although the high-redshift observations remain very pre- 
liminary, the low-redshift results seem to be converging 
toward a low escape fraction of ionizing radiation, less than 
10%, not only from galaxies such as the Milky Way but 
from starburst galaxies. The discrepancy between the 
observations and the theoretical predictions of starburst 
galaxies may be because the gasdynamics driven by star- 
bursts are not considered in the theoretical models. In this 
paper, we study the impact of supershells and galactic winds 
on the escape of ionizing radiation from dwarf starburst 
galaxies. We focus on dwarf galaxies because they are 



relatively simple systems. In addition, small-scale structures 
such as dwarf galaxies are the first ones to collapse in the 
early universe, according to the theory of hierarchical 
galaxy formation, and therefore may make a significant 
contribution to the reionization of the universe. 

To study the escape fraction of ionizing radiation from 
dwarf starburst galaxies, we model disks with baryonic 
masses Mj = W, lO^, and Mg at redshifts z = 0, 3, 
and 8. Our motivation is to study the impact of supershells 
and galactic winds due to starbursts on the escape of ioniz- 
ing radiation by using the Mac Low & Ferrara model (1999; 
hereafter MF99), combined with the photoionization code 
of Abel, Norman, & Madau (1999). We assume disks that 
are rotationally supported with exponential surface density 
profiles, on the basis of the predictions of standard hierarch- 
ical cosmologies (Mo, Mao, & White 1998). We note that 
the presence of a well-formed disk is uncertain at high red- 
shift (e.g., Abel et al. 1998; Abel, Bryan, & Norman 2000) 
but argue that such idealized disks yield a lower limit to the 
escape fraction of ionizing radiation from a halo with given 
baryonic mass. In addition, our model is limited to a single 
ionizing source (i.e., starburst clump or star cluster) at the 
center of the disk, but we argue that centrally located 
sources also yield a lower limit to the escape fraction. In § 2, 
we describe our photoionization model with Lyman 
continuum photon luminosities, based on the Starburst99 
model (Leitherer et al. 1999). In § 3, we describe our disk 
and halo models and how we model star formation in the 
disks. In § 4, we summarize our numerical methods, devel- 
oped to model the impact of repeated supernovae (SNe) in 
the galactic disks, based on MF99. We show results on the 
escape fractions from undisturbed disks without supershells 
in § 5.1 and those from disks disturbed by starburst bubbles 
with supershells and outflows in § 5.2. In § 6, we compare 
our model results with observations of dwarf starburst gal- 
axies and also attempt to compute the contribution of dwarf 
starburst galaxies to the UV background radiation at high 
redshift. Conclusions foUow in § 7. 



2. PHOTOIONIZATION MODEL 

We compute the propagation of ionization fronts by using 
a photon-conserving radiative transfer code, developed to 
simulate inhomogenous reionization self-consistently within 
a cosmological hydrodynamical simulation (Abel et al. 
1999). It was slightly modified for our problem. The code 
solves for the time at which the ionization front (I-front) 
reaches each grid cell along many different rays by solving 
the Stromgren equation along each ray: 

dri 

47rri — whi = N^h{t) - 47raB j ripHer^ dr , (1) 

where rj is the radial distance of the I-front at that angle 

from an ionizing source, «h i is the number density of neu- 
tral hydrogen, rip and «g are number densities of protons 
and electrons, Nph{t) is a photon luminosity, and is the 
case B recombination constant. The on-the-spot approxi- 
mation is used to treat the effects of the diffuse emission of 
ionizing photons. This is expected to be reliable for 
scenarios in which stellar sources are responsible for the 
ionization of hydrogen (Abel et al. 1999). 

The code solves for the radiative transfer around a point 
source; in our study this is a central starburst source. We 



52 



FUJITA ET AL. 



Vol. 599 



choose a single site of star fomiation for rmmerical simplic- 
ity but also based on the assumption that some sort of 
dynamical instability, accompanied by the transfer of 
angular momentum, drives gas to the center. However, local 
observations of dwarf galaxies show that the starbursts 
occur in isolated clumps in the disks where the column 
densities of gas are greater than 10^' cm"^ (Skillman 1997), 
so our assumption of a central source can certainly be 
improved. Another way of looking at the problem is that we 
do not actually calculate the escape fraction of ionizing 
radiation from a realistic starburst galaxy but rather from a 
single star-forming clump or star cluster in the galaxy. We 
later argue in § 3 that we are computing the lower limit for 
the escape fraction for the single clump by placing it at the 
center and that the result with a point source can be 
extended to a case in which there are multiple ionizing 
sources scattered around the disk. 

We consider two different starburst scenarios: (1) an 
instantaneous starburst and (2) multiple starbursts spread 
over a time of 20 Myr. The photon luminosity below the 
Lyman limit of the cluster is taken from the Starburst 99 
model (Leitherer et al. 1999), as a function of time. The 
model is based on a power-law initial mass function (IMF) 
with exponent a = 2.35 between low-mass and high-mass 
cutoff masses of A/|ow = 1 M--., and Mup =100 Mr..,, respec- 
tively, with metallicities ranging from iF = 0.001 to 0.04. 
This approximates the classical Salpeter (1955) IMF. Most 
observations of star-forming and starburst regions are con- 
sistent with the Salpeter IMF, although uncertainties can be 
large (Leitherer et al. 1999). In Figure 1, we plot the photon 
luminosity as a function of time A^ph(f) for three different 
models in which 10^ of stars between 1 and 100 Mq 
form, producing ~40,000 SNe. Note that with the instant- 
aneous-burst model, most photons are produced in the first 
few megayears of the burst, during the lifetime of the most 
massive stars. For different strengths of starburst, we line- 
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Fig. 1. — ^Number of Lyman continuxim photons produced per second 
as a function of time -/Vph(?). Three different models are shown: the 
instantaneous-burst model (solid line), multiple-burst model convolved 
with a timescale of 20 Myr (dashed line), and constant star formation model 
with a duration of 100 Myr (triple-dot-dashed line). 



arly scale the above luminosity, A'^ph(/, M = 10^ Mq), 
according to the mass available for the starburst, MgB, as 
follows: iVph(/, Msb) = iVph(f, 10^ M^M^/i^ x 10^ Mg), 
with the coefficient, ^ = 2, which accounts for the mass of 
stars below 1 A/;, with the power-law turnover below 
M = 0.1 Mq. The number of photons produced per solar 
mass is 5.8 x 10^". 

The radiation transfer code is a postprocessing step that 
operates on a given density distribution at a given time step 
in our simulation. Therefore, the time evolution of bubble 
structures (see § 4) does not include ionization. However, 
the additional pressure due to photoionization is small com- 
pared with the bubble internal pressure. We can accurately 
compute the instantaneous escape fraction /esc(0 (see § 5.1 
for the definition) at a given time, since the ionization 
front propagates sufficiently rapidly to adjust almost 
instantaneously to a changing density distribution. 



3. DISK 

The escape fraction is very sensitive to the H i distribution 
in a galactic disk. In our study, we assume that there is a 
single ionizing source at the center of the galactic disk. In 
this case, the shortest axisymmetric path for the photons to 
escape is in the vertical direction. With the same amount of 
gas and a fixed column density distribution from the mid- 
plane, Nii,{R) = Jq°° nui{R, Z)dZ in axisymmetric cylin- 
drical coordinates, it is harder for ionizing photons to 
escape in the vertical direction from a thick disk with a large 
scale height H than from a thin disk with a small H. This is 
because the gas located high above the midplane at distance 
Z can trap the photons reaching that height, which are 
diluted by the inverse square law, 1 /AttZ^, more efficiently 
than the gas located close to the midplane, which encounters 
a larger photon flux, despite the higher density closer to the 
midplane. We show later in this section that the formation 
of shells of swept-up ISM by superbubbles complicates this 
situation. 

To show that it is harder for photons to escape a thick 
disk, it is useful to define the effective recombination rate in 
a given direction through a disk: 



a-^HpHer 



(2) 



where r = {R?- -I- Z^)'^^ is the distance from the central ion- 
izing source located at f = 0. The photon luminosity per 
steradian along each ray at a given time, N^\,{t)/A'K, must 
exceed H for any photon to escape in the given direction. 

For example, if we assume the vertical distribution of 
neutral hydrogen gas is Gaussian with number density 
«(Z) = exp(— or exponential with «(Z) = 
exp(— Z///), the radial column density from the ionizing 
source at R = Q is iVHi,o = A'hK^ = 0) « (V7r/2)«*//' 
(Gaussian) or n*/f (exponential), so the effective 
recombination rate in the Z direction is 



1 

8V2 
1 



for Gaussian 
for exponential 



(3) 
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When A^Hi,o is fixed, a larger scale height H means a larger 
effective recombination rate Sz. 

On the other hand, if we assume that the radial distribu- 
tion of gas is exponential with number density n{R) = 
n^. exp{—R/ Rd), where Rj is the scale radius for column 
density distribution, then the effective recombination rate in 
the R direction is 



-2 n3 



(4) 



When A'hi.o is fixed, a large H means a smaller effective 
recombination rate, H«. It is easier for the photons to move 
in the radial direction in a thick disk with a larger scale 
height, H, than in a thin disk with a smaller H. Note that 
Er > Ez as long as Rd > H, although the radial density 
profile in our model disk falls off faster than exponential. In 
the next subsection, we constrain Rd and Nwi{R) by a 
standard theory of hierarchical galaxy formation. 



3.1. Disk and Halo Model 

We set up a disk with an exponential surface density 
profile. Eg = So exp(— where Rd is a scale radius, 
consistent with observations of local compact dwarf gal- 
axies (van Zee, Skillman, & Salzer 1998). To model a disk 
with an exponential surface density profile, we follow Mo 
et al. (1998), who compute the scale radius Rd by assuming 
that a fraction of halo mass, = Md/M/,, settles into a 
non-self-gravitating, rotationally supported disk in a singu- 
lar isothermal halo with a fraction of halo angular momen- 
tum, j'^ = Jd/Jh- The halo angular momentum is produced 
by tidal torques from the background density distribution. 
The conventional way to specify the angular momentum is 
by the dimensionless spin parameter, \ = Jh\E\MJ^^^^, 
where J/, is the total angular momentum of the halo and E is 
its binding energy. A^-body experiments show that the distri- 
bution of A is Gaussian, with a peak at A ~ 0.05 (e.g., Cole 
& Lacey 1996). 

If gas with given angular momentum settles into an 
exponential disk, the scale radius is given by 



Rd = -^ ( — |Ar200 



(5) 



where r2oo is the virial radius of a halo, defined as the radius 
at which the mean overdensity in the halo is 200pbg and p^g 
is the background density. The central surface density is 
then given by Eq = Md/lnRj, which enables us to fully 
specify the shape of the disk. This method provides a good 
fit to the observed size distribution of galactic disks with 
Jd/md = 1 (Mo 1998). 

Equation (5) is derived with an isothermal profile for the 
halo, 



Phir) 



(6) 



despite the unphysical singularity at its center, since with 
this profile the total binding energy E can easily be calcu- 
lated analytically. The central structure does not markedly 
influence Rd- Once we have defined Rd using this approxi- 
mation, though, we use a softened isothermal halo with a 
central core to specify other halo parameters and to solve 
for the galactic structure in hydrostatic equilibrium. We 
use the form given by Salucci & Persic (1997) based on 



observations of dwarf galaxies, 
Phir) 



Po 



1 + (r/ror 



(7) 



where ro and po are the core radius and the central density 
(Binney & Tremaine 1987). The virial radius is then 



1/2 



(8) 



As in MF99, we calculate ro and po as a function of halo 
mass based on the empirical relation between ro and po for 
nearby dwarf galaxies (Burkert 1995) extrapolated to the 
early universe by using the usual cosmological scaling 
relations. 



ro = (8.9 X 10"* kpc) 



1 Mp 



1/2 



-1/3 



PO = (6.3 X 10'" Mq kpc"') 



1 Ma 



-1/3 



no 



(9) 



(10) 



We extended the relations to high-redshift galaxies by 
assuming that the ratio of the core radius to the virial radius 
and the ratio of the core density to the background density 
stay constant. 

The disks in our models are set up in hydrostatic equili- 
brium in the potentials of the dark matter halos and, in 
some models, the disks themselves. The gas that is impor- 
tant for trapping ionizing radiation is all assumed to be in 
the warm neutral state, heated and stirred up by star forma- 
tion prior to the major starburst that we model, yielding an 
effective sound speed, including a turbulent contribution, of 
Cj ~ 10 km s"' with corresponding effective temperature 
T 10^-^ K. The situation is more complicated in a real 
multiphase gas disk. However, only a few SNe from massive 
stars suffice to blow through a thin disk of cold neutral 
medium similar to those observed (e.g., Mac Low & 
McCray 1988), thereby creating a funnel for ionizing radia- 
tion to escape. In addition, ionizing radiation escapes 
through a thin disk with a small scale height more easily 
than a thick disk with a large scale height. Therefore, we do 
not consider the thin cold disk in our calculations. 

MF99 used the dark-to-baryonic mass ratio cp = Mh/ Mg 
observed in local spiral galaxies, 4> — 34.7Af~7'^^, where 
Mgj = Mg/10'^ Mq„ hereafter where "Icl" indicates 
"local" (Persic, Salucci, & Stel 1996). The relation was 
derived from a sample of larger galaxies with Mh > 10'" Mq 
in the local universe, but MF99 extrapolated the result to 
lower masses because the trend of increasing dark matter 
fraction with decreasing galaxy mass is indeed consistent 
with some observations (e.g., Mateo 1997). However, the 
extrapolation of this result to higher redshift is imperfect. 

In Figure 2, we plot Mi, as a function of Mg expected with 
the observed relation and with the standard picture of 
the hierarchical ACDM cosmology by using JIq = 0.37 and 
= 0.05, which corresponds to = 7.4. The difference 
appears in both lower and higher mass disks. Using the (/)ici 



54 



FUJITA ET AL. 



Vol. 599 



10'2 



10" 



r 






Persic et ol. (1996) : 




ACDM model 







10' 108 109 IQIO 10" 

log (Mbisk/Msun) 



Flo. 2. — Mass of dark matter halos as a function of the mass of baryonic 
matter in the halos, with the relation (Persic, Salucci, & Stel 1996) and 
^ = lA expected with ACDM model (Hq = 0.37 and fij = 0.05), where 
4> = Mu/Mj. 

relation underpredicts Mg in lower mass halos and over- 
predicts it in higher mass halos, compared with the value of 
Mg predicted in the ACDM model, in which baryons are 
assumed to simply follow the dark matter. The under- 
prediction can be interpreted to mean that a significant 
amount of gas was depleted or blown away in the shallow 
potentials of small dwarf galaxies because of the past star- 
burst activities, as seen in the simulations of MF99. The 
overprediction at high masses can be interpreted to mean 
that in larger galaxies the deeper potentials retain more gas 
in the halos and the cooling is also more efficient (ocw^). 
However, we show in § 5.1.5 that the difference in Mg 
between (/)ici and (p = 7.4 is too small to influence the escape 
fractions of ionizing photons in any of the halos we 
examine. 

The scale heights of the disks are determined by setting 
the gas into hydrostatic equilibrium in the halo and, in some 
cases, the disk potentials. When a disk potential is included, 
we use the thin-disk approximation (Toomre 1963). The 
additional disk potential makes a disk thinner than other- 
wise. We refer to the disks with both dark matter and disk 
potentials as " thin " disks and those with only dark matter 
potentials as " thick " disks. We list the scale heights of our 
model disks in Table 1. The vertical gas distributions are 
roughly exponential in the thin disks but Gaussian in the 
thick disks. The presence of disk potentials makes the disks 
unreasonably thin when large surface densities are expected 
at high redshift and in large galaxies. For example, a thin 
disk with M^ = Mq at z = 8 has a scale height ^ of less 

TABLE 1 

Scale Heights of Thin and Thick Disks 



Mj ThinDisk,z = Thick Disk, z = z = 3 z = 8 

(Mq) (pc) (pc) (pc) (pc) 



10^ 100 580 71 21 

10' 56 750 93 27 

IQi" 19 980 121 36 



than 10 pc. Therefore, we consider our thick disks to be 
more appropriate to describe high-redshift disks, especially 
when highly dynamical processes of cooling and heating are 
expected after halo formation. A thin disk with Mj = 10'" 
Mq at z = has only H = 19 pc, but a thick disk has 
H = 980 pc, which is unreasonably large. Therefore, we 
experiment with both thick and thin disks as the opposite 
extrema of possible gas distributions within the observa- 
tional and theoretical limits. 

We do not consider preexisting dumpiness of the disk, 
but we do consider the extreme case of dumpiness: shell for- 
mation due to starburst bubbles. If shells of swept-up ISM 
are formed at r^^^ with average density «sh ^nd thickness Ar^h, 
the column density from the central ionizing source is 
A'sh = WshA/'sh in spherical coordinates [r = {Z'^ + R^Y^'^] 
and the effective recombination rate is 

Ssh = OEnshA^sh^'sh , (11) 

according to equation (2). Since the preferred direction for 
the growth of bubbles and therefore for the escape of pho- 
tons is the vertical direction, we compare the effective 
recombination rate due to the shell, E^h oc UshNshZj^^, with 
that in an undisturbed disk, Ez oc N^^qH, from equation 
(3) in the vertical direction. We assume a Gaussian or 
exponential distribution of hydrogen as before. 

We can assume that the shell, at least initially, is bounded 
by an isothermal shock because the swept-up gas in the 
shells cools very effectively in our model galaxies, as we 
show in § 5.2. In the absence of magnetic fields, the isother- 
mal shock compression is nsh = J^^n{Zsh)- A perpendicular 
magnetic field, however, reduces the compression behind 
an isothermal shock to «sh = \/2-^A«(-Zsh), where .Mi^ = 
!;(47rp)'/^B~' is the Alfvenic Mach number of the shock. 
In the case of a Gaussian atmosphere, H(Zsh) = 
«* exp(-Z2j^///2), while in an exponential atmosphere 
«(Zsh) = «* exp(— Zsh///), where «* is the central number 
density. The ratio of the effective recombination rate due to 
the shell to that in an undisturbed disk is 

„ ^^„^ _ 1 2\/2«sh(Zsh) for Gaussian 
\ 8nsh(Zsh) for exponential 

in the limit of Zjh > H. The shells trap the ionizing photons 

more effectively than the undisturbed disk, as long as Z^\,IH 
is less than [log(2y2,y//^)] ''^^ or log(8,/#^) , respectively, or in 
the magnetized case (log4.y//A)'/^ or log {%\fi.Mp^, respec- 
tively. For example, is greater than Hz if the bubble is 
within 2.4/f in a thick disk {H = 580 pc) and within 7.0/? in 
a thin disk (/f = 100 pc) with M = 10^ Mq at z = 0, with 
.M = 10. 

In addition, the starburst-driven shells can trap the radia- 
tion more effectively than numerous dumps distributed 
throughout the disk, because such clumps would leave only 
a diluted ISM and allow for a larger escape fraction (Wood 
&Loeb 2000). 

Our standard model applies a cold dark matter (CDM) 
cosmology with = 1 and h = 0.5 and the observed dark- 
to-baryon ratio ^i^, as in MF99, and spin parameter 
A = 0.05. We experiment with both thin and thick disks, 
that is, both including and neglecting disk potentials. We 
also show in § 5 . 1 . 5 that changing the cosmology to the stan- 
dard ACDM (f2o = 0.37, = 0.63, h = 0.7, = 0.05) 
does not influence the escape fractions. 
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There are limitations to our study. We model idealized 
rotationally supported disks with uniform single-phase ISM 
and exponential surface density profiles. Our model disks 
approximate the dwarf galactic disks observed in the local 
universe reasonably well, but we do not know whether such 
disks initially form at high redshift after halos collapse. A 
reahstic scenario is that the stars form as the gas cools and 
collapses then explode as SNe and stir up the rest of the 
cooled gas in the halos. Disks may not form at all before star 
formation can occur, according to the study of the forma- 
tion of first stars with high-resolution three-dimensional 
adaptive mesh refinement simulations by Abel et al. (1998, 
2000). We think that a more turbulent dynamical ISM than 
the one we study will allow more photons to escape through 
funnels of hot gas. In this study, we solve for the lower limits 
of escape fractions by examining well-formed rotationally 
supported disks at all redshifts. 

3.2. Modeling the Star Formation Rate 

We considered several scenarios for the star formation in 
our disks since few empirical constraints are available for 
high-redshift galaxies. These scenarios share some common 
characteristics. First, we use a threshold density to delimit 
the star-forming region of the disk. This constraint effec- 
tively differentiates the star formation in systems with 
different angular momentum. Second, we assume that the 
amount of star formation is proportional to the amount of 
gas in this region. And, finally, merger-driven star forma- 
tion is recognized by considering gas distributions that are 
more centrally concentrated than our fiducial disk model. 
We define a star formation (or starburst) efiiciency param- 
eter and define the amount of gas available for a starburst 
asMsB =hMd. 

The density threshold for star formation in nearby spiral 
galaxies (Kennicutt 1989) and starburst nuclei (Kennicutt 
1998) is well described by models for the onset of gravita- 
tional instability (Kennicutt 1989, 1998). The critical surface 
density for star formation is determined by a column density 
cutoff based on Toomre's criterion, E^. = O.SkCs/ttG, where 
K is the epicyclic frequency. The coefficient 0.5, hereafter the 
gravitational threshold, is fitted to the observations by 
Martin & Kennicutt (2001). We assumed velocity disper- 
sions maintained by turbulence, Cj = 10 km s~', to get the 
highest efficiency in trapping radiation and so to yield a 
lower limit to the escape fractions. This assumption is con- 
sistent with the values of velocity dispersions measured in 
nearby galaxies (van der Kruit & Shostak 1984; Dickey, 
Hanson, & Helou 1990; van Zee et al. 1998) but is uncertain 
at high redshift. It is also not clear that the same relation 
holds in dwarf galaxies (Hunter, Elmegreen, & Baker 1998). 
To explore other cutoffs, we also considered a constant 
column density cutoff, = 10^' cm"^ (Skillman 1997). 

The critical masses of our disks with Md =10^ Mg at 
z = 0, 3, and 8 are/*Mrf with/* = 0.0056, 0.097, and 0.32, 
respectively, by the gravitational threshold and /* = 0.013, 
0.79, and 0.94, respectively, by the 10^' cm^^ cutoff. Our 
standard model for the gas distribution produces fairly low 
star formation rates in the disk with Md =10^ Mg at z = 0. 
The observations of blue compact galaxies (BCDs) by van 
Zee et al. (1998) show that there is another population of 
dwarf galaxies with higher central column densities and 
these tend to be the bursting systems most likely to contrib- 
ute to escaping UV photons. Therefore, we increased the 



CN 
1 


23 




Rd=0.78 kpc at z=0 


E 






■ 


o 


22 






X 

z 






: 
- 




21 






cn 
O 




"■:r7^V^.^;^j7ii£;z^- - — . _^ 


CM 

1 






— \ — , — , — , — , — \ — , — , — , — , — 
Rd=0.28 kpc at z=3 ; 


E 








o 


22 






X 

z 








o 


21 






O 




CM 

1 


23 




Rd=0.12 kpc at z=8 


^cm" 


22 










column density ! 


X 

Z 






\ gravitational threshold - 


o 

CT> 

O 


21 




102icm-2 cut-off 



0.0 0.5 1.0 1.5 

R (kDc) 



Fig. 3. — Surface density profiles for the disks with = Vfi Mg at 
z = 0, 3, and 8, compared with the critical surface densities predicted by the 
gravitational threshold and the 10^' cm^^ cutofT. The sizes of the disks are 
with /Jrf oc (1 + z)"' . The critical surface density by the gravitational 
threshold is not drawn at the center because the gas is supported by pressure 
rather than rotation there. 



central gas density in our model disk with Md = 10^ Mq at 
z = by decreasing the spin parameter A to 0.035 (within 
1 a deviation of the mean, A = 0.05). The new critical mass 
is f^Md with /* = 0.074 with the gravitational threshold 
and/* = 0.2 with the 10^' cm"^ cutoff. Figure 3 shows the 
surface density distributions of the disks with Mj = 10** Mq 
at z = 0, 3, and 8 with the critical surface densities predicted 
by the gravitational threshold and the 10^' cm"^ cutoff. In 
this study, we consider three possible star formation efficien- 
cies for single starburst clusters in our model disks: 
/* = 0.006, 0.06, and 0.6. Note that /* also depends on the 
amount of metals present in dwarf galaxies, but the effects 
of metallicity will be smaller than the 2 orders-of-magnitude 
difference in/* . 

Figure 4 shows the distributions of exponential scale radii 
Rd and central surface densities Eq for our model galaxies 
based on the standard model, compared with those for five 
BCDs by van Zee et al. (1998) and a number of late-type 
dwarf galaxies by Swaters et al. (2002). The scale radii and 
the surface densities of our model galaxies with Md =10* 
Mq and A = 0.035 and with M,, = 10' A/., and A = 0.05 at 
collapse redshift z = are consistent with the observed dis- 
tributions. We cannot in principle compare our higher red- 
shift galaxies with the local BCDs and late-type dwarf 
galaxies. We also note that the CDM predictions of the disk 
spin parameter are very uncertain (Dalcanton, Spergel, & 
Summers 1997). A disk with a low spin parameter has 
a higher surface density and so a higher effective 
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Fig. 4. — Distributions of exponential scale lengths Rj and central sur- 
face densities So in units of H i column density (cm^^) of our model galaxies 
at z = (triangles), r = 3 (squares), and z = 8 (diamonds) and those based 
on the observations of blue compact dwarf galaxies (BCDs; Fig. 1 1 from 
van Zee, Skillman, & Salzer 1998) and of late-type dwarf galaxies (Swaters 
etal. 2002). 

recombination rate: Ez oc A^^^q oc A""*. However, the same 
disk yields a higher critical mass and so a higher star forma- 
tion efficiency/* oc A^^'* based on the Schmidt law (Schmidt 
1959): M* oc Eq The increase in Ez is compensated for by 
the number of ionizing photons A^ph due to the increase in 
/* . Therefore, the quahtative results on the escape fractions 
of ionizing photons in a disk with a given baryonic mass can 
be applied to the same disk with a different spin parameter 
A. Despite the uncertainty in A distribution, our models at 
z = trace the upper boundary of the properties of local 
dwarf galaxies well (Fig. 4). The actual form of high-redshift 
disks remains the largest uncertainty, as we mentioned in 
the end of § 3.1. 

A last limitation of our model is that we deal with only a 
single star-forming source at the center of a disk. In realistic 
starburst galaxies, we expect that multiple starburst clumps 
are scattered around the disks (e.g., Vacca 1996; Martin 
1998). Recall that the photons find their way out of the disks 
first in the vertical direction (Ez < Eji) and that 
Ez oc N^X-^)H{R) and A^hi(^) oc exp{-R/Rd), although 
H increases as a function of R. This means that a clump 
placed at the center yields the lower limit for the escape frac- 
tion and the escape fraction we compute is a lower limit for 
the same starburst clump placed off center. A starburst 
clump forming above Z — has a higher escape fraction on 
the side it forms. The lower limits of escape fractions 
from reahstic galaxies can be inferred by summing the 
contributions from starburst clumps in the disks. 



4. SHELL FORMATION AND GALACTIC OUTFLOWS 

Large outflows of gas are observed from both dwarf and 

massive starburst galaxies in the local universe (e.g. Heckman 
et al. 2000; Martin 1999; Lehnert & Heckman 1996) and at 
high redshift (Franx et al. 1997; Pettini et al. 1998). These 
outflows are driven by the mechanical energy supplied by 
massive stars in the form of SNe and stellar winds (Tomisaka 



& Ikeuchi 1988; De Young & Heckman 1994; Suchkov et al. 
1994, 1996; MF99; D'Ercole & Brighenti 1999; SiUch & 
Tenorio-Tagle 1998, 2001). Deep Ha imaging often shows 
filamentary or bubble-like structures in dwarf galaxies 
(Marlowe et al. 1995; Delia Ceca et al. 1996; Martin 1998), 
which are thin dense shells of swept-up ISM. When there is 
sufficient time for the shells to cool and become dense, such 
cooled neutral shells of hydrogen may be able to trap the 
ionizing radiation inside the host halos. 

Dove et al. (2000) studied the escape of ionizing radiation 
from the disk of the Milky Way and found that the shell of 
the superbubble driven by each OB association quickly 
traps or attenuates the ionizing flux from that association. 
The escape fraction of photons does not increase even after 
superbubbles of larger associations blow out of the H i disk 
and form dynamic chimneys because by that time the pro- 
duction of Lyman continuum photons has drastically fallen. 
However, the escape fraction was already very small, ~0.1, 
in their Galactic disk without these shells. This is because 
they modeled the disk as having OB associations with a 
power-law photon luminosity function, and the largest OB 
association has a photon luminosity of only 2 x lO'i s~^ 
Even in dwarf starburst galaxies, much higher photon lumi- 
nosities are expected, greater than a few times 10^^ s~', 
although they probably come from multiple star clusters 
(e.g., Martin 1998; see § 6.1). The collective power of stellar 
winds and tens of thousands of SNe in starburst galaxies 
can be much larger than that from any single OB association 
studied in their model. In our study, we examine how super- 
shells and galactic outflows influence the escape of ionizing 
radiation from powerful starburst clumps. 

To model the effects of shell formation due to repeated 
SN explosions from dwarf starburst galaxies, we use ZEUS- 
3D, a second-order Eulerian astrophysical gasdynamics 
code (Stone & Norman 1992; Clarke 1994) that uses 
Van Leer (1977) monotonic advection. Runs were done on 
Silicon Graphics Origin 2000 machines using eight process- 
ors and typically took ~ 10-1 5 hr in the models at z = for 
the computation up to 10 Myr and about 4 days in the 
models at z = 8 for computation up to 5 Myr after the onset 
of starbursts. 

Our numerical methods and initial conditions generally 
follow MF99 (see their § 4): we use axisymmetric ratioed 
grids; a cooling curve by MacDonald & Bailey (1981), line- 
arly scaled by a factor of 10 to account for low metallicity in 
dwarf galaxies; and a tracer field (Yabe & Xiao 1993) to 
avoid overcooling in the bubble interiors. Note that the cool- 
ing curve A(r) in principle does not linearly scale; however, 
we also tested the cooling curve of Sutherland & Dopita 
(1993) with primordial gas composition. It agrees with the 
linearly scaled A{T) by MacDonald & Bailey (1981) within a 
factor of a few in the temperature range between 10^ and 10^ 
K, which are the shock-heated temperatures of the shells 
expected in our simulations with Vsh « 100 km s~^ We also 
find that the cooling of sheUs A{T)n^ is not so sensitive to 
A{T) but primarily to the density of the shells, n. 

However, there are a few changes. As in MF99, the 
thermal energy source that drives a constant luminosity 
wind is set up as a source region with a radius of five zones. 
Since the sizes of galaxies decrease as (1 + z)~\ the resolu- 
tions adopted also decrease roughly as (l+z)~'. This 
means that physical sizes of the source regions decrease at 
high redshift. We justify this by noting that smaller source 
regions favor shell formation. Here we test the maximum 
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effect of shell formation on the escape of ionizing radiation, 
and smaller source regions give a firm lower limit. 

The source region is set up initially in pressure equili- 
brium with the background, but at a temperature a factor r] 
larger than that of the background and with density corre- 
spondingly lower by a factor of I/77. The source region then 
has energy added to it at a rate E = Lmech/ V and mass at a 
rate M = poE/{rieo), where po and eo are the midplane mass 
and energy density of the background disk, respectively, to 
roughly keep the specific energy of the source region con- 
stant. The amount of mass in this wind accounts not only 
for SN ejecta but also mass that would have evaporated 
from the shells if we had included thermal conduction in our 
model. In fact, evaporation from the shells provides most of 
the mass within a realistic bubble (Weaver et al. 1977). 
MF99 fixed r] = 1000, which was appropriate for the small 
range of mechanical luminosities Zmech that they studied. 
However, to take into account our broader range of L^sxh, 
we vary r] as a function of Lmech by directly equating the 
mass input rate M with the mass evaporation rate predicted 
by Weaver et al. (1997) for a spherical bubble. 



Me 



,27/35 4/35 -2/35 



(13) 



As we have t and Lmech, we can estimate rj by noting the 
weak dependence of the mass evaporation rate on the den- 
sity and choosing to set the density constant hq ~ 10^ cm^^. 
In this way, we keep the density and temperature within the 
bubble fairly close to the analytic solution of Weaver et al. 
(1977) and so can follow the interior dynamics with some 
degree of rehability. 

We note that photoevaporation by stellar UV radiation 
from interior molecular clouds increases the amount of 10"* 
K gas within a superbubble. This will increase the surface 
area available for conductive flow, although it is hard to 
estimate how much will end up in the hot interior gas 
(McKee, van Buren, & Lazareff 1984; ShuU & Saken 1995). 
However, within the range of mechanical luminosities con- 
sidered in our study, Lmech ^ 10^^ ergs s"', the cooling due to 
the photoablated material of ~ 10"* Mq will not be important 
after the bubble grows greater than ~30 pc because of 
reduced interior density (ShuU & Saken 1995). 

We list the resolutions used for the computation of the 
escape fractions of ionizing radiation in undisturbed disks 
and in disks disturbed by galactic outflows in Table 2. We 
find that we need to resolve the scale heights of the undis- 

TABLE 2 

Minimum Grid Resolutions Used in Modeling 
THE Disks 



Mj z = z=3 z = 8 

(Mq) (pc) (pc) (pc) 

10« 2.5 0.63 0.28 

10' 4.4 1.1 0.49 

10'» 7.5 1.9 0.77 

Note. — The miniirmm grid resolutions used in 
modeling the disks are roughly proportional to 
Ryir <x mI'^{1 + z)~\ We selected all the disks at 
z = and the disks with Mj = 10* Mq at z = 3 
and 8 for the bubble experiments (first row and 
first coliunn). We also experimented with higher 
resolutions, 1.25 and 0.38 pc, for the model with 
Mj = 10* Mq at z = to test the sensitivity of 
shell cooling to the resolution. 



turbed disks with more than 50 zones to get converged 
results with our photoionization code. For the outflow 

models, we are not able to fully resolve the shell formation 
when the shocks are isothermal with shell densities 
«sh = itd^'^ or «sh = nj\/2./M in the presence of magnetic 
fields, where na is the density of the ambient disk. 

The cooling time of the shells formed by initially adia- 
batic shocks is tc = 3kTs/nA{T), with postshock electron 
number density He = n. The shocked temperature is 



7; = i3fijl6k)vl ~ 1.3 X 10= 



Vsh 

100 km s- 



K, 



(14) 



where = (14/23)oth is the mass per particle in the fully 
ionized postshock gas, assuming «He/«H = 0.1, and fsh is 
the shock velocity. To estimate the cooling time we use 
the analytic estimate of the Gaetz & Salpeter (1983) cooling 
function between 10^ and 10^ K from Mac Low & McCray 
(1988): A(r) = (1.0 x 10-23 ergs cm^ s-^)CT^'^^ with = 
T/IO^ K, where C is the metallicity relative to solar ( = 
^/^Q. Then we find that the instantaneous cooling time of 
the postshock gas in the shells is 



tc = 1.2 X lO"* 



yr 



100 km s- 



3.4 



1 cm" 



O.I 



(15) 



which remains short (tc < 1 Myr), while the bubble is within 
a few scale heights of the disk, with the expected Vsh^200 
km s Therefore, we expect isothermal shocks near the 
disks. 

We show in Figure 5 an example of our resolution study 
of the behavior of isothermal shocks, for a thin disk with 



108 Mp, at 



0, with a mechanical luminosity 

10"*" ergs s"'. The expected density of swept-up 
behind unmagnetized isothermal shocks is 



z = 

-1 



-^mech ■ 

shells 

«sh = Ud^^. The effective sound speed = 10 km s"' in our 
simulations is driven by turbulence and photoionization, so 
the Mach number could be larger given the short timescales 
on which turbulence dissipates and hot stars disappear. 
Magnetic fields on the order of a microgauss may, on the 
other hand, reduce the compression in realistic galaxies to 
«sh = ndy/2Ji\, as much as an order of magnitude lower. 




Fig. 5. — Resolution study for the disk with Af^ = 10* Mq at z = with 
mechanical luminosity Lmech = 10*° ergs s"'. The density of the shells is 
extremely sensitive to the resolution of the simulations. The theoretical 
value rid Jl^ is indicated. 
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Figure 5 shows that the unmagnetized isothermal shell 
structure is better resolved as the resolution is increased 
from 25 to 0.38 pc, but the expected peak density of the 
shells due to isothermal shocks remains about a factor of 4 
larger than the peaks in our simulations. At the same time, 
the fragmentation of the shells due to Rayleigh-Taylor 
instability is better resolved in the simulations with a higher 
resolution. The shells created by isothermal shocks would 
trap the ionizing radiation more effectively (S a «sh), but 
the fragmentation of the shells opens up holes through 
which the ionizing radiation can escape. 

Note that the development of Rayleigh-Taylor instability 
in our model is limited not only by the resolution but also by 
the assumption of azimuthal symmetry. Mac Low, McCray, 
& Norman (1989) argued that the assumed symmetries of 
the problem largely determine the specific modes of the 
instability that are excited and thus the distribution of mass 
into particular fragments (e.g., rings in azimuthal symme- 
try) but that the appearance of unstable regions is independ- 
ent of the symmetry assumed. We expect more and smaller 
shell fragments in three dimensions, but the escape of ioniz- 
ing radiation is primarily determined by the opening angles 
of the chimneys created after bubbles blow out, as we show 
in § 5.2, so the suppressed modes of shell fragmentation 
within those angles will not substantially affect the outcome 
of our computations. We also do not expect thermal insta- 
bility (Chevalier & Imamura 1982) or linear dynamical thin- 
shell instability (Vishniac 1983) to be important for shell 
fragmentation, based on the numerical simulations of these 
instabilities by Blondin & Cioffi (1989) and Mac Low & 
Norman (1993), respectively. 

To estimate the actual ability of the shells to trap ionizing 
radiation in our Umited resolution simulations, we set the 
shell density to the analytic value for an isothermal shock 
and compute the recombination measure Hsh at the posi- 
tions of the shells that are not fragmented because of the 
instability. (This correction is needed only for disks at z = 
because even unresolved shells in high-redshift disks are 
dense enough to trap the ionizing photons.) To make the 
correction, we set a temperature threshold Tth = 5 x lO"* K 
in our photoionization code, to identify cooled shells that 
have not yet fragmented. Then we make a correction 
Ssh = a^naJ^^Nshrl^^ according to equation (11) if the maxi- 
mum density in the identified shell is smaller than 
Hd-Jf^ ~ nd[vsh/{lO km s~^)]^, the approximation applied 
within the uncertainties of magnetic fields and thermal 
properties of the shells. This approximation tends to overes- 
timate the shell density at larger scale heights in low-redshift 
galaxies but is appropriate to test the maximum effect of the 
shells on the escape of ionizing radiation. We can reliably 
read A'sh from the lower resolution simulations, since it is 
independent of cooling. Note that shell formation is more 
sensitive to the resolution than to the cooling function A{T) 
applied with linear metallicity dependence because the 
cooling is proportional to \{T)rP-. 

We assume mechanical luminosity to be constant in time 
for a given /* in our simulations, by taking the average of the 
luminosity =Sfmech(0 drawn from the Starburst 99 model 
(Leitherer et al. 1999). We apply J" = 0.008 at z = based 
on the observations of local dwarf galaxies (Skillman, 
Kennicutt, & Hodge 1989; Hunter, Gallagher, & 
Rautenkranz 1982) but = 0.001 at high redshift, assuming 
that the timescale for metal replenishment by previous star 
formation prior to the major starbursts that we model is very 
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Fig. 6. — Mechanical luminosities -2'mccii(0 is a function of time, drawn 
from the Starburst 99 model (Leitherer et al. 1999) when 10^ M^, of gas is 
converted to stars, for the instantaneous-burst model with metallicities 
Z = 0.008 (thick solid line) and = 0.001 {thin solid line) and for the 
multiple-burst model with = 0.008 (thick dashed line) and = 0.001 
(thin dashed line). 



short, much less than the Hubble time at a given redshift 
(Wada & Venkatesan 2003; A. Fujita et al. 2003). Figure 6 
shows the evolution of .^mech(f) as a function of time when 
10^ Mq of gas is converted to stars based on four different 
starburst scenarios. When the metallicity is fixed, the 
mechanical luminosity i^mech(0 is about an order of magni- 
tude lower in the first few megayears with the multiple-burst 
model than with the instantaneous-burst model, although 
^mech(0 soon riscs to the same average mechanical luminos- 
ity expected with both star formation models. When the star- 
burst model is fixed, ifmech(0 is about an order of magnitude 
lower in the first few megayears with iF = 0.001 than with 
iF = 0.008, during the time S'mech{t) is dominated by stellar 
winds. This occurs because stellar wind strength depends 
strongly on metallicity. The mechanical luminosity =Sfniech(^) 
with 3^ = 0.001 soon approaches ^mech(0 expected with 
iF = 0.008. In this study, we are specifically interested in the 
evolution of bubbles at early time before they blow out of the 
disks, at / < 10 Myr at z = and t < 3 Myr at z = 3 and 8 (see 
§ 5.2 for discussion on bubble dynamics). 

The average mechanical luminosity Z-mech within 10 Myr 
after the onset of starburst with Msb = 10^ M,-., is Lmech(10^ 
Mq) f» 2 X 10"*'' ergs s"' for the instantaneous-burst model 
with iF = 0.008 (see Fig. 6, thick solid line). This is used as a 
reference value for the other three cases. We lower Lmech by 
1 order of magnitude from the reference value if the multi- 
ple-burst model is applied with the same metallicity and if 
the instantaneous-burst model is applied with iF = 0.001. 
We lower Lmech by 2 orders of magnitude from the reference 
value if the multiple-burst model is applied with ^ = 0.001. 
Note the weak dependence of bubble radius on mechanical 
luminosity in a uniform medium; 7?bi oc /-^jf^h (Castor, 
McCray, & Weaver 1975). We vary the mechanical luminos- 
ity linearly with the total mass of gas converted to stars, 
Msb: Z.mech(MsB) = L^^^{IQ^ Mq)Msb/(^ x 10^ Mq) with 
e = 2. 
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We model the effects of starbursts on dwarf galactic disks 
with Md = 10^ 10^ and lO^o at z = and with 
Md =10^ Mr., at z = 3 and 8. Because of limited computa- 
tional resources, we simulate the evolution of bubbles only 
until they blow out of the disks and create galactic outflows. 
Our aim is (1) to test the effects of supershells on the entrap- 
ment of photons while the shells are within the disks, that is, 
before the blowouts, and (2) to test whether galactic out- 
flows increase the escape of photons after the blowouts. The 
opening angles of chimneys created by the outflows are 
assumed to stay constant after they blow out. 



5. RESULTS 

5.1. Undisturbed Disks 

We first analyze the results for disks without outflows. 
We define the total escape fraction of photons from a disk 
as 



(16) 



where A^ph(0 is a given photon limiinosity, with the 
instantaneous escape fraction 



/esc(0 



1 



/7r/2 /•27r 
dsined6 / d(p 
■jr/2 Jo 



aBn^ir, 6, ip)r^dr 



(17) 



The last integration gives the total effective recombination 



rate, where ri is the radius of the ionization front at time t. 
Neither the total nor the instantaneous escape fractions 
depend on the viewing angle of the observer. It is important 
to remember that observers measure the instantaneous 
escape fraction fesM, not the total escape fraction F^.^^. 
We show later in this section that /esc(0 and F^sc can differ 
drastically when the Lyman continuum spectrum depends 
strongly on time. 

Figure 7 shows the fraction of photons escaping from 
disks as a function of star formation efficiency /* with our 
standard model. With the maximum star formation effi- 
ciency that we consider, = 0.6, which roughly corre- 
sponds to the critical mass set by the 10^' cm^^ cutoff, most 
of the photons ( > 50%) can escape in all the disks at all red- 
shifts with any given spectrum of Lyman continuum radia- 
tion, except in thick disks at z = 8. The total escape 
fractions /esc decrease as a decreasing function of /* but vary 
with the different assumptions for disk morphology and 
star formation history that we apply. We discuss below 
the dependence of various parameters in our study on the 
escape of ionizing photons. 

5.LL Scale Heights 

The thick disks can trap the ionizing radiation more effec- 
tively in many cases, especially at high redshift, while thin 
disks let more radiation out. This is because the photons 
find their way out of the disks in the vertical direction first, 
while they are often trapped in the radial direction 
(Hz < Er). The effective recombination rate in the vertical 
direction is Ez oc H, from equation (3). Recall that the pres- 
ence of disk potentials makes the disks unreasonably thin at 
high redshift; therefore, the results with thin disks should be 




Fig. 7. — Total escape fractions Fesc for models without bubbles with thin (filled symbols) and thick (open symbols) disks and with the instantaneous-burst 
model (diamonds and solid lines) and the multiple-burst model (triangles and dashed lines). The total escape fractions decrease as redshift increases and do 
not depend much on the star formation model applied. Thick disks can trap more ionizing photons than thin disks in general. 
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ignored at z > but can still be compared with the results 
from thick disks. There are some cases in which the thick 
disks yield slightly higher /esc than the thin disks, when the 
ionization front immediately propagates through the entire 
disks, while Er > Ez and Er oc H-^ from equation (4). 

5.1.2. Redshift and Mass 

The escape fractions decrease as redshift increases in 
undisturbed disks (see Fig. 7), confirming previous work by 
Wood & Loeb (2000) and Ricotti & ShuU (2000). The trend 
is due to the increase of the surface density, T, oc R^^ oc 
(1 + zf. Recall Sz oc R^'^H oc (1 + z)^H, so the effect of an 
increase in the redshift is substantially greater than that of a 
decrease in H (see Table 1). However, we show in § 5.2 that 
the redshift dependence for a given f^f no longer exists once 
the effects of galactic outflows are taken into account. On 
the other hand, we do not see the trend of decreasing escape 
fractions as the disk masses increase, as in Wood & Loeb 
(2000). The escape fractions at a given redshift are approxi- 
mately the same for all our disks with a given because 
Ez oc Eg// oc M2/3//andAfph oc/*M. 

5.1.3. Star Formation History 

There is a mild increase in the escape fractions from the 
multiple-burst model to the instantaneous-burst model. The 

trend reflects the behaviors of photon luminosities N^Y^t), 
which are an order of magnitude higher during the first few 
megayears with the instantaneous-burst model than with 
the multiple-burst model (see Fig. 1). We show in Figure 8 
the instantaneous escape fraction /esc(0 and the cumulative 
total escape fraction Fg^^J^t) as a function of time in a disk 
with Md = 10>* Afg at z = with/* = 0.006 and 0.06 and 
two different models of Lyman continuum history. The 
shape of /esc(0 as a function of time reflects that of A'^ph(0- 
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Fig. 8. — Instantaneous escape fractions yisc(0 (thin lines) and the total 
cumulative escape fractions F^(t) (thicic lines) plotted for a model with 
=10* Mq at z = with Lyman continuum histories, based on the 
instantaneous-burst model and the multiple-burst model. The star 
formation efficiencies /* applied are 0.06 (solid and dashed lines) and 0.006 
{triple-dot-dashed line and dotted line). The figure shows that the observed 
value of the escape fraction at time t,f^{t), does not correspond to the total 
escape fraction F^sc = lim,^oo Pex{t). 



For example,/.sc(f) peaks sharply in the instantaneous-burst 
model within the first few megayears, as A'^pii(0 does, and as 
a result most Lyman continuum photons escape the disk 
quickly, yielding the total escape fraction /"esc ^0.8. On the 
other hand, /sc(0 rises slowly to its peak in the multiple- 
burst model, yielding F^sa « /esc(0 < 0.8. We emphasize 
that the instantaneous escape fraction /esc(0 that observers 
can measure does not necessarily reflect the total fraction of 
ionizing photons that escape out of starbursting galaxies 
/'esc if the bursts are nearly instantaneous. 

/'esc and/.sc(f) with a different star formation history can 
be inferred based on the results above. The peak photon 
luminosity primarily determines the escape fraction from a 
starburst clump. Therefore, F^^^. and fesdt) from multiple 
starbursts with a longer duration or a constantly star- 
forming clump with a correspondingly lower peak photon 
luminosity are expected to be smaller. However, they do not 
scale linearly with the peak photon limiinosity, because of 
disk geometry. 

5. 1 .4. Distribution of Ionizing Sources 

Recall that our study is limited to single starburst clusters 
or ionizing sources at the center of our model disks. These 
central clusters convert /*Afrf of gas into stars, and so /* in 
our study does not correspond to star formation efficiencies 
of the entire disks. We also argued in § 3.2 that we compute 
lower limits to the escape fractions in our study by placing 
the starburst clusters at the center and that the escape 
fractions remain as lower limits when there are multiple 
starburst clusters of the same luminosities scattered around 
the disks. Below, we give some insights into a situation in 
which a given amount of gas is available for star formation 
throughout a given disk. 

Dove & Shull (1994) and Dove et al. (2000) modeled the 
Galactic disk as having randomly distributed OB associa- 
tions. The total escape fraction from the Galactic disk is 
found to be -0. 1 . In the study by Wood & Loeb (2000), star- 
burst ionizing sources are spread throughout the disks with 
emissivity scaling as (like OB associations) or «hi, and 
quasar sources are located at the center of the disks. At 
z > 5, more than 0.3 of Lyman continuum photons escape 
from their miniquasars and none escape from their starburst 
galaxies. The low escape fractions seen in some of their 
model disks are due to the distribution of ionizing sources. 

In a given disk, there is a minimum photon limiinosity 
^ph,min for any photon to escape or for a Stromgren sphere 
to break out of it, as is also mentioned in Dove & Shull 
(1994). This means when a fixed amount of mass is con- 
verted to stars in a galaxy, having numerous or multiple ion- 
izing sources is not a favorable condition for the photons to 
escape, since many OB associations then have photon Ivmii- 
nosities below A'ph.mm- In Table 3, we hst A'^phmin for our 
model disks. Note that A^ph.mm is equal to the minimum 
effective recombination rate, that is, Ez oc N^^H, times 47r. 
Therefore, #ph,miii is higher for more massive disks and disks 
at high redshift and is also higher for thicker disks. We con- 
clude that when a fixed amount of mass is available for star- 
burst, placing it all in a single source gives the most 
favorable condition for the escape of ionizing radiation. 

5.1.5. Otiier Issues 

Changing the cosmology from an SCDM model (f2o = 1 
and h = 0.5) to a ACDM model {Uq = 0.37 and h = 0.7) 
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TABLE 3 

Minimum Photon Luminosity for a Photon to 
Escape from Model Disk 
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Note. — The minimum photon luminosity for any photon to escape from 
the model disk is -/Vph,mm = 47rHj. 



changes the virial radius of the halos, 7?vir oc 
[h'^VLfi/Q,{t)Y^^^ , and therefore the scale radius Rd oc R^r 
and the central surface densities of the disks Eq oc R'^^. At 
high redshift f2(/) ~ 1, so the changes in and fio cancel 
out. Any significant dilference between cosmologies appears 
only for low-redshift disks at 1 with Vl{t) = fio- Even 
there, the maximum increase in Eq is only a factor of 2. We 
find in our low-redshift ACDM models that the qualitative 
results on the total escape fractions do not change. 

We also checked the effect of increasing the disk gas in 
low-mass halos at all redshifts by using the cosmological 
dark-to-baryonic ratio (j) = 1 A instead of A high value 
of means a decrease in halo mass for a given disk mass. The 
virial radius and the scale radius of the new halo decrease 
{Rd oc and the surface density of the disk increases 

up to a factor of a few. We find that there is no significant 
change in /'esc from our model disks with 4> = 1 A. 

The surface density increases only by a factor of a few 
when the spin parameter of a halo differs by 1 cr from the 
mean: Eq oc A-^ = (0.035/0.05)"^ « 2. The above results 
with different cosmology and values of (j) suggest that the 
change in the angular momentum of a disk within a. \ a 
range will not influence the qualitative trends of /"esc seen in 
Figure 7. Any change in A much greater than 1 a should be 
compensated for by the change in star formation efficiency 
/c , as we discussed in § 3.2. 

If we use the dark matter profile by Navarro, Frenk, & 
White (1997) with the halo concentration factor, c = 5-30 
(Mo 1998), instead of the dark matter profile in equation 
(6), the maximum decrease in scale radius Rd is expected to 
be less than j. The increase in Ez is then limited to less than 
1 order of magnitude. 



5.1.6. Summary 

We can draw our results from Figure 7 together by con- 
sidering the effective recombination rate in the vertical 
direction from equation (3): Hz oc Rd'^H oc (1 + z)^//. The 
total escape fraction of ionizing radiation /"esc decreases as 
redshift z increases. For a disk with given mass Md and 
column density distribution E(7?) at a given redshift, F^^c 
also decreases as the scale height H of the disk increases. 
Changes in scale radius Rd due to changes in cosmology, 
dark-to-baryonic matter ratio, spin parameter A by ~ 1 ct, or 
dark matter distribution produce only small changes in 
-Fesc- As we argued in § 3.2, a further increase (decrease) in 
Rd, for example, with a greater than 1 a increase (decrease) 
in A should be compensated with a decrease (increase) in/^it 
and A'^ph. Therefore, we conclude that the escape of ionizing 
radiation is sensitive to vertical gas structures rather than 



radial gas structures for fixed Md, specifically the scale 
height H. 

Despite the sensitivity of the escape fraction to disk mor- 
phology, the results converge at low redshift. At z = 0, most 
of the photons easily escape even with star formation effi- 
ciency/* as low as ~0.006. Recall we are solving for the 
lower limits on the escape fractions for starburst clusters. 
Our results strongly suggest that most photons do escape 
from dwarf galactic disks at low redshift if bubble dynamics 
are not considered. These results seem to contradict the 
observations of local dwarf starburst galaxies with very low 
escape fractions, /esc(0 < 0.06. However, we show in the 
next section that the formation of shells by starburst 
bubbles significantly alters the results. 

5.2. Disks with Shells and Galactic Outflows 

We model the effects of stellar winds and SNe from star- 
bursts in our dwarf galactic disks to examine how swept-up 
shells and galactic outflows affect the escape of ionizing 
radiation. We experiment with disks with Md = 
108-lOiOMg at z = and Md = 10^ Mg at z = 3 and 8. We 
show below two significant results of our experiments: (1) 
swept-up shells of neutral hydrogen can substantially trap 
the ionizing radiation before blowouts, even in disks at low 
z that otherwise yield high /esc(?) and F^sc and (2) galactic 
outflows allow a significant fraction of ionizing radiation to 
escape directly to the IGM from disks at high redshift that 
otherwise yield Fg^c ~ 0. 

5.2.1. Low-Redshift Thin Disks 

Figure 9 shows the propagation of ionization fronts over- 
laid on density distributions in a thin disk with Md =10^ 
Mq at z = 0, with = 0.06 by using the multiple-burst 
model (with ^ = 0.008) at ? = 2 and 6 Myr after the onset 
of the starburst. We take -Lmech = lO"*" ergs s^', 1 order of 
magnitude lower than the mechanical luminosity expected 
with the instantaneous-burst model (see Fig. 6), to examine 
the effects of shells swept up by a bubble before it begins to 
blow out (/<5Myr). 

The swept-up shell has density nd{R,Z)\/2./M \ < 
^sh < nd{R,Z).yM^ because of the isothermal shock. As we 
mentioned in § 4, we compensate for the inability to fully 
resolve the density jump in the shells in our simulations by 




Fig. 9. — Entrapment and the escape of ionizing radiation shown in black 
on top of density distributions in a thin disk with = 10* Afg at z = 0, at 
t = 2 and 6 Myr after the onset of starburst with /* = 0.06. The Lyman 
continmun history applied is the multiple-burst model with S = 0.008, 
which corresponds to an average mechanical luminosity L^sA = 10** ergs 
s~' within 5-10 Myr. The cold dense shells of swept-up ISM completely 
trap the ionizing photons at ? = 2 Myr {left) but begin to fragment as the 
bubble accelerates by f = 6 Myr (right). The fragmentation creates funnels 
for the photons to escape to the IGM. 
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Flo. 10. — Changes in the instantaneous escape fractions fox(t) when the 
effect of shells of swept-up ISM is included. The model is a disk with 
=10* Mq at z = 0, with /* = 0.06, and the Lyman continuum photon 
luminosities are based on the instantaneous-burst model {diamonds and 
solid line) with L^ech = 10*' ergs s~' and the multiple-burst model {triangles 
and dashed line) with Lmech = 10'"' ergs s~'. The instantaneous escape frac- 
tionsy^sc(') with shells cooled to the density expected for isothermal shocks 
are shown with open symbols, /esc(0 with unresolved shells in our lower 
resolution simulations are shown with filled symbols, andyisc(') without 
shells are shown with Unes. The decrease in /esc(0 is significant once the 
maximum effect of fully cooled shells is included, before the bubble blows 
out of the disk at « cs 8 Myr. 

computing the eifective recombination rate, Ssh oc «siiA'sii''sh' 
at each angle through the shells in our photoionization 
code, using the hydrodynamical value for «sii to maintain a 
lower limit. At ? = 2 Myr in Figure 9 (left), the dense shell 
completely traps the ionizing radiation, so /esc(0 = 0. By 
t = 6 Myr (right), the bubble starts to accelerate in the verti- 
cal direction, where the disk gas is most stratified, so the 
Rayleigh-Taylor instability begins to fragment the shells. 
The fragmentation opens funnels for the photons to escape, 
so /esc(0 = 0.31. By that time, the accelerated shock in the 
vertical direction is no longer isothermal. The further evolu- 
tion of the bubble leads only to an increase in fesc{t), because 
the bubble eventually blows out with more fragmentation 
and the mechanical luminosity of the burst is expected to 
increase up to lO'*' ergs s~' after f ~ 10 Myr. 

Figure 10 shows the instantaneous escape fraction fesdt) 
for the model shown in Figure 9 with = 0.06 with shells 
cooled to the theoretical density «sh = n^.^^ (open symbols), 
with the unresolved shells from our simulations (filled sym- 
bols), and without shells (lines). We compute _/esc(0 every 
2 Myr up to 10 Myr in our simulation with a resolution of 
2.5 pc. The photon luminosity of the instantaneous-burst 
model reaches its peak, ~2 x 10^^ s^', already at t = 0.2- 
0.3 Myr and decreases by more than an order of magnitude 
by ; ~ 5 Myr. That of the multiple-burst model is ~2 x 10^^ 
s~^ by f ~ 2 Myr, still rising very slowly up to ~3 x 10^^ s~' 
for the duration of starburst (tsb = 20 Myr). 

A significant effect from the underdensity of the shell due 
to numerical resolution appears only at t<A Myr, before 
the bubble begins to blow out of the disk. The shells can trap 
the photons very effectively while the bubble is within a few 



scale heights of the disk (H = 100 pc). Afterward the shells 
fragment as the bubble accelerates, allowing the photons to 
escape. Figure 10 shows that the proper treatment of shell 
formation is crucial to study the full effect of the shells on 
the escape of ionizing photons. 

With the fully cooled shells, we significantly suppress 
/esc(f)- As a result, the total escape fractions F^sc decrease 
from 0.95 to 0.16 with the instantaneous-burst model and 
from 0.78 to 0.42 with the multiple-burst model. We assume 
that /esc(0 stays constant after Z > 10 Myr because the open- 
ing angle for the escape of photons is primarily determined 
by that of the bubble after blowout, which remains fairly 
constant at ~30°^0°. The value of Fgsc computed with the 
multiple-burst model is strictly a lower limit, since the 
mechanical output of the burst increases by an order of 
magnitude at later time (see Fig. 6). However, recall that the 
bubble evolution depends only weakly on Lmecii: ^bt 
-^mech- decrease in Fesc is more significant with the 
instantaneous-burst model than the multiple-burst model. 
This is because the shells significantly trap ionizing photons 
only while the bubble is still within a few scale heights, and 
most photons are produced during that time with the 
instantaneous-burst model. 

5.2.2. Low-Redshift Thick Disks 

Bubble evolution depends sensitively on the scale height 
Hin our model disks. The results for /esc(0 and are dif- 
ferent in a thick disk with a larger H. As mentioned above, 
the shells of swept-up ISM are very effective in trapping 
radiation before acceleration starts above a few scale 
heights. The time when the bubble reaches H in the vertical 
direction can be approximated as ?h « ^mlWo^^^^^ ^ 
uniform medium with a constant density, mq (Castor et al. 
1975). We can approximate hq = «* ~ ^Hi q/H, that is, the 
central density of the Gaussian or the exponential profile. 
Then for constant column density A^hi,o, tn oc jj^g 
scale height of the thin disk with dark matter and disk 
potentials is 100 pc, but that of the thick disk with only dark 
matter potential is 580 pc, both with Md = 10^ at z = 0. 
Therefore, the time for the bubble to reach // is a factor of 
10 shorter in the thin disk than in the thick disk. In addition, 
the vertical gas distribution is exponential in the thin disk 
because of the additional gravitational pull by the disk, 
while it is Gaussian in the thick disk. The stratification 
drives the acceleration of a bubble (Kompaneets 1960; 
Mac Low & McCray 1988), so the bubble blows out of the 
thin disk much earlier than the thick disk. On the other 
hand, the shock velocity t;sh at H is proportional to 
'^mecii-'^ '''^' '■^^ Mach number is a factor of 2 larger in 
the thin disk than in the thick disk. Recall that Ssh oc Hsh 
and «sii = nj./M- , but this decrease in Ssh is small compared 
with the increase in the time for a blowout in the thick disk. 
The opening angles of the chimneys cleared out by 
Rayleigh-Taylor instability is ~30°^0° at the time of 
blowout in both disks. 

Figure 1 1 shows that the instantaneous escape fractions 
_/esc(f) are significantly suppressed in the thick disk with a 
larger scale height, until the bubble begins to accelerate at 
Z « 10 Myr. The instantaneous escape fractions /esc (0 are 
sensitive to the different scale heights, H, because blowouts 
determine the transition from the entrapment of photons by 
shells to the escape of photons through chimneys. The com- 
puted total escape fractions F^^ are 0.16 in the thin disk and 
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out. The scale height of the disk strongly determines the 
time of blowout. The scale heights of the thin disks are 56 
and 19 pc and those of the thick disks are 750 and 980 pc 
with Mj = 10' and 10'" M^„ respectively. Therefore, the 
bubbles in thick disks can trap the radiation up to ^ ~ 10 
Myr after the onset of starbursts, but the bubbles in thin 
disks immediately blow out of the disks and create chimneys 
for the escape of photons. As a result, the total escape frac- 
tions F^sc are significantly more suppressed in thick disks 
with large H than in the thin disks with small H (see Table 4) 

5.2.3. High-Redshift Disks 

At high redshift, halos are very dense [p/, oc (1 + zf] and 
disks are compact [Rj oc (f) ' and Eq oc (1 +2)^]; there- 
fore, the strong gravitational field of the halos creates highly 
stratified atmosphere in the disks even without disk poten- 
tials (thick disks in our study). As a result, bubbles quickly 
blow out of the disks at very early time, about a few mega- 
years after the onset of starbursts, although the mechanical 
power of the bursts is expected to be smaller, while stellar 
winds dominate in high-redshift disks with ^ « 0.001 than 
in low-redshift disks with ^ « 0.008. 

Figure 12 shows the evolution of a starburst bubble and 
the propagation of ionization fronts at ? = 2 and 3 Myr in a 
disk with Md = 10^ at z = 8, with = 0.06 and the 
multiple-burst model. We expect lower metallicity at high 
redshift, .^<0.001; therefore, the star formation efficiency 
/* = 0.06 with the multiple-burst model corresponds to 
J^mech = 10^' ergs"' withiu the first few Myr. This is because 
the stellar winds dominate the mechanical luminosity dur- 
ing the time, and the production of stellar winds is very sen- 
sitive to the metallicity (see Fig. 6). We apply Lmech = 10^' 
ergs s~' because we are specifically interested in the 
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Fig. 11. — Instantaneous escape fractions /isc(0 in the thick disk with 
H = 580 pc (bold diamonds and bold triangles'), compared with those in the 
thin disk with if = 100 pc (diamonds and triangles) for the same galaxy 
shown in Fig. 10. The thick disk traps more photons than the thin disk until 
i 10 Myr, because the growth of the bubble in the thick disk is slower 
than that in the thin disk. 



0.11 in the thick disk with the instantaneous-burst model 
and 0.42 and 0.12, respectively, with the multiple-burst 

model. 

Experiments with high-mass disks with Md =10' and 
10'° at z = with = 0.06 show that the qualitative 
results are the same: /esc(?) is suppressed only before blow- 



NoTE. — Total escape fractions F^^ are computed in the disks with supershells and galactic outflows created by 
starburst bubbles (shell), compared with F^^ in undisturbed disks (none). Only thick disks are considered at z > 0. The metal- 
licities assumed are Z = 0.008 at z = and ^ = 0.001 at z > 0. 
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Fig. 12. — Starburst bubble in a disk with M = 10* MQ at ; = 8, with = 0.06. the multiple-burst model, and a corresponding mechanical luminosity of 
lO'"' ergs s"'. The bubble sweeps the high-redshift dense ISM into very dense shells up to r = 2 Myr (left), which are very elfective in trapping the ionizing 
radiation. The bubble, however, immediately blows out of the disk hy t = 3 Myr (right). The chimneys created by the outflow let ~20% of ionizing photons 
produced by stars escape to the IGM. 



evolution of bubbles at early time before they blow out of 
the disks and the shells fragment. Recall that we experiment 
with thick disks only at high redshift. Figure 12 shows that 
the bubble immediately breaks out of the disk by ? = 3 Myr, 
even with the very small mechanical power supplied by low- 
metallicity stellar winds, and lets the photons escape 
through the chimney created, so /esc(0 = 0.014 at ? = 2 
Myr, but fesc{t) = 0.24 at / = 3 Myr. The opening angles of 
the chimneys cleared out by Rayleigh-Taylor instability are 
'-~^40°-50° at Z = 3-5 Myr. The fragmented shells are already 
falling back and prevent the escape of photons; however, 
after t = 3 Myr, more and more supernovae go off, supply- 
ing more mechanical power for the bubble to grow (see Fig. 
6), so we expect the escape fraction to remain /esc(z) <;0.2. 
Remember that only a negligible fraction of photons can 
escape out of undisturbed high-redshift disks because of 
their high density. 

Figure 13 shows the increase in fesc{t) when the effects of 
galactic outflows are taken into account in disks with 
Md = 10^ Mq at z = 3 and at z = 8 with /* = 0.06 and 
Lmtch = 10^" ergs s"' for the instantaneous-burst model and 
^mech = lO-'^ ergs s"' for the multiple-burst model. The peak 
luminosities are the same as before. At z = 8, the escape 
fractions fescit) are nearly zero with both Lyman continuum 
histories in undisturbed disks, but/esc(f) increases to greater 
than 0.2 as the bubbles blow out of the disk at a very early 
time, by ? ~ 1 Myr (instantaneous) and ~3 Myr (multiple). 
At z = 3,/esc(f) are rather high with the instantaneous-burst 
model but are less than 0.15 with the multiple-burst model 
in undisturbed disks, but /esc(0 are greater than 0.2 as the 
bubbles blow out by Z ~ 3 and ~6 Myr, respectively. We did 
not further compute the evolution of the bubbles after they 
blow out, because of the great computational expense. 
However, we do not need to compute fescit) after blowouts 
if the starbursts are instantaneous. When the starbursts last 
for 20 Myr, we assumed that the sizes of the chimneys for 
the escape of photons do not change much after blowouts. 

We find that early blowouts in dense, high-redshift disks 
create chimneys for the photons to escape freely to the 



IGM. As a result, the total escape fractions /'esc are 0.28 
(0.056 in the same but undisturbed disk) with the instanta- 
neous-burst model and 0.22 (0.0) with the multiple-burst 
model in the disk at z = 8, and are 0.13 (0.44) and 0.20 
(0.11), respectively, in the disk at z = 3. The total escape 
fraction Fg^c is suppressed with the instantaneous-burst 
model at z = 3 because most ionizing photons are produced 
before the blowout, which occurs slightly later at z = 3 than 
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Fig. 13. — Instantaneous escape fractions foso(t) in undisturbed disks 
(line) and in disks with the effects of shells and galactic outflows (symbols) 
as a function of time. The models are disks with Mj = 10* M0 at z = 3 
(bottom) and z = 8 (top), with /* = 0.06. The Lyman continuum histories 
are based on the instantaneous-burst model, which yields imech = 10"" ergs 
s~' (solid line and diamonds), and on the multiple-burst model, which yields 
imech = 10^' ergs s~' (dashed line and triangles). The escape fractions/esc (?) 
are computed up to when the bubbles blow out of the disks and galactic 
outflows are formed. Note that /esc(?) increases once the outflows create 
funnels for the ionizing photons to escape directly to the IGM (except in the 
disk at z = 3 with the instantaneous-burst model). 
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at z = 8. We conclude that the dynamic chunneys created 
by galactic outflows increase the escape of ionizing radia- 
tion in the high-redshift disks, in which otherwise only a 
negligible amount of photons can escape. With Fesc^O.2, 
dwarf galaxies may make a substantial contribution to 
reionizing the IGM in the high-redshift universe. 

5.2.4. Other Issues 

We discussed only the results with f^^ = 0.06. However, 
we can guess the trend of/esc(0 and T^esc as a function of/>it, 
since the formation of dense swept-up shells due to isother- 
mal shocks that trap ionizing radiation is expected for all 
the disks with mechanical luminosities corresponding to 

< 0.6 at all redshifts, as long as the bubbles are within a 
few scale heights of the disks. The star formation efficiency 

has a linear dependence on Z^mech, so the time for a blow- 
out, /bi oc ^~ech' changes by only a factor of 2 with an 
order-of-magnitude increase in/* and L^nech- In addition, 
the photon luminosity Np\i{t) has a linear dependence on 
while the effective recombination rate of the shells, 
Ssh oc n^hM^, is only a factor of ~5 larger with an order- 
of-magnitude increase in A'ph(?), and Lmech because 
fsh oc Therefore, an early increase in /esc(0 and an 

increase in Fesc are expected with/* = 0.6. The opposite can 
be said of the behaviors of /esc(?) and F,.^ with/* = 0.006. 
The summary of the effects of shells and galactic outflows 
on the total escape fractions Fesc with different /* is listed in 
Table 4. 

We did not run the simulations for large disks at high red- 
shift, because the extremely strong shocks in high-density 
environments are unstable in ZEUS-3D. However, the 
extrapolation of our results indicates that the blowouts 
occur very quickly in high-redshift massive disks, so the 
total escape fraction is increased to .fesc 5:0.2. Note that 
the redshift dependence of Fesc discussed in § 5.1 does not 
exist in the presence of galactic outffows. 

We do not consider the effects of dust on the escape frac- 
tions of photons in our models. However, we think that the 
qualitative results will not change in our models even if real- 
istic distributions of dust are included because the dust con- 
tent is expected to be low in low-metallicity systems such as 
dwarf galaxies, especially at high redshift, and the dust will 
also be swept away with the galactic outflows. Note also 



that a hydrogen column density of 10^' cm-^ is required for 
dust to produce raust = 1 even with solar metallicity (Spitzer 
1978, eq. [7-23]), while an H i column of only 1.5 x 10'^ 
cm^^ produces thi = 1 at the Lyman edge. Thus, dust can 
be neglected unless hydrogen is extremely ionized. 

6. DISCUSSION 
6. 1 . Comparison with Observations 

In this study, we examined the effects of supershells and 
galactic outflows on the escape fraction from a single star- 
burst cluster at the center of our model disk. Although this 
is only a first approximation to a realistic dwarf starburst 
galaxy, we can gain some insight by comparing our results 
with the observations of three dwarf starburst galaxies with 
escape fractions less than 0.06: NGC 1705, NGC 5253, and 
NGC 4214. In particular, observers measure the instantane- 
ous escape fraction fesdt), not the total escape fraction Fgsc- 
The value of /esc(0 depends on starburst age and the stage 
of bubble evolution, as shown in the previous section, and 
therefore may not be representative of /'esc- 

First, we go over the basic properties of the three dwarf 
starburst galaxies. Table 5 gives a summary. 

1. NGC 1705 is a "nucleated" blue compact dwarf gal- 
axy with multiple ellipsoidal shells expanding at Vsh ~ 100 
km s~^ driven by a super-star cluster NGC 1705-1 (Meurer, 
Freeman, & Dopita 1992; Meurer et al. 1995; Marlowe, 
Meurer, & Heckman 1999; Heckman et al. 2001). The ioniz- 
ing continuum of the order of 10^^ s"' is produced mainly 
by the background star-forming population, not by the 
nucleus, although it dominates the UV-optical flux from the 
galaxy. The nucleus has mass -^1.5 x 10^ with an esti- 
mated age of ~13 Myr after a nearly instantaneous burst, 
and star formation rate in the background is 0.03 Mq yr~' 
with a star formation timescale of ~1 Gyr (Meurer et al. 
1992). The H i rotating disk shows a hint of disturbance 
with the optical/H i peaks offset from the H i dynamical 
center, but it is not clear whether the starburst was exter- 
nally triggered (Meurer, Staveley-Smith, & Killeen 1998). 

2. NGC 5253 is an amorphous dwarf galaxy with two 
Wolf-Rayet clusters. These clusters produce most of the 
observed ionizing radiation ~4 x 10^^ s~' (Schaerer et al. 



TABLE 5 

OnsLRVhi) PKOFhRiiLs oi NGC 1705. NGC 5253. AM) NGC 4214 



Parameter NGC 1705'' NGC 5253'^ NC;C4214'' 



rf(Mpc) 6.1 4.1 3.6 

Mb (mag) -16.20 -17.62 -17.65 

H I extent (kpc) ^-4.8 4.3x4.3 14.6x12.9 

Pciic(kms~') -60 38-60 70 

«(fcirc) (kpc) 2.1 1.24 2.91 

Type Amorphous/BCD Im(pec), amorphous/BCD IAB(s)m 

Miyn(MQ) 1-2E9 >4.3E8 3.3E9 

Mh,(M0) 8.7E7 1.4E8 1.8E9 

e(s-i) -1E52 3.86E52 1.94E52 

/esc(%) <5.8 <6.0 <4.0 



Note. — The escape fractions y^sc(0 ii these dwarf galaxies are estimated to be less than 0.06 
(Heckman etal. 2001). 

" The observed properties of NGC 1705 are taken from Meurer et al. 1992, 1998 and Marlowe 
etal. 1999. 

*" The observed properties of NGC 5253 and NGC 4214 are taken from Tables 1, 4, and 5 in 
Martin 1998. 
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1997) and drive superbubbles observed in Ha probing fila- 
mentary structures expanding at ~35 km s, coincident with 
soft X-ray emission (Martin & Kennicutt 1995; Martin 
1998; Marlowe et al. 1999). The clusters are very young, 
with estimated ages 2.8 and 4.3 Myr after (nearly) instanta- 
neous bursts and are embedded in an older star-forming 
background (0.1-1 Gyr). The measured stellar rotation rate 
of only less than 7 km s"' around the optical minor axis and 
the high stellar velocity dispersion (50 km s"') suggest that 
the galaxy was a dwarf elliptical before the current starburst 
began (Caldwell & Phillips 1989). An H i study of NGC 
5253 shows that the bulk of the neutral hydrogen gas rotates 
about its optical major axis, and the peculiar H i kinematics 
may be a signature of recent accretion, disruption of tidal 
companions, or interaction with the nearby spiral galaxy 
M83 (Kobulnicky & SkiUman 1995). 

3. NGC 4214 is a (barred) irregular dwarf galaxy, also 
listed as a Wolf-Rayet galaxy. The total ionizing continuum 
of 3 X 10^^ s~' is mainly produced by several very young 
(~3 Myr) massive star clusters, which also drive superbub- 
bles sweeping the ISM into shells with > 50 km s~', 
creating two large H ii complexes (Martin 1998; MacKenty 
et al. 2000). It has a very extended H i disk with a peak 
column density ~3 x 10^^ cm~^, and the overall H i distri- 
bution indicates the existence of H i holes and shells, also 
with a hint of a bar (Mclntyre 1997; Walter et al. 2001). The 
total H I mass and estimated dynamical mass are larger than 
those of NGC 1705 and NGC 5253 (Martin 1998; see 
Table 5). 

In NGC 5253 and NGC 4214, the ages of the super- 
clusters are estimated to be <4 Myr. As shown in Figure 11, 
the instantaneous escape fraction/esc (0 is suppressed to less 
than ~0.1 before the bubbles start to accelerate at ? « 4 
Myr in the thin disk or at ? « 10 Myr in the thick disk with 
Md =10* Mq at z = 0. Our low-redshift disk models pre- 
dict low escape fractions in NGC 5253 and NGC 4214 
because of their young age. The observed kinematics sup- 
port the hypothesis that the shells have not broken up yet. 
The low escape fraction measured in NGC 1705 requires a 
different explanation since the age of the supercluster is ~ 1 3 
Myr. However, multiple shells extending out to ~1 kpc 
from the supercluster do not yet seem to be fragmenting 
because of Rayleigh-Taylor instability (Meurer et al. 1992). 
The low escape fraction in NGC 1705 is partly explained by 
our results that the blowout time is delayed until ~10 Myr 
in thick disks with large scale heights. 

Escape fractions have been observed for other galaxies: 
blue compact galaxies with Mb > -20 such as Mrk 66, Mrk 
1267 (Leitherer et al. 1995; Hurwitz et al. 1997), and Mrk 54 
(Deharveng et al. 2001), nucleated starburst galaxies, IRAS 
08339+6517 and Mrk 496 (Leitherer et al. 1995; Hurwitz 
et al. 1997), Lyman break galaxies (Steidel et al. 2001), and 
high-redshift galaxies from the Hubble Deep Field (HDF) 
with magnitudes Mi5ooa~ —21 to —25 (Fernandez-Soto 
et al. 2003). The escape fractions from the low-redshift 
BCDs and starburst galaxies are estimated to be less than 
10%. On the other hand, the escape fractions from LBGs at 
z ~ 3.4 from Steidel et al. (2001) are estimated to be greater 
than ~50%, and those from HDF galaxies at 1.9 < z < 3.4 
from Fernandez-Soto et al. (2003) are estimated to be less 
than ~15%. Although most of them are larger galaxies in 
which starburst activities extend to rather more than a few 
kiloparsecs (e.g., Giavahsco, Steidel, & Macchetto 1996), a 



simple extrapolation of our results with dwarf galaxies sug- 
gests that the escape fractions should correlate with the 
stages of bubble evolution and therefore with the ages of 
starbursts in the galaxies. 

For the low-redshift BCDs and starburst galaxies, the 
ages of starbursts are estimated to be a few to ~6 Myr, and 
these starburst clusters seem to drive large-scale motions of 
the interstellar gas at fsh about several hundred km s^^ 
(Gonzalez Delgado et al. 1998). In light of our results, the 
swept-up shells around the young starburst bubbles may 
trap the ionizing photons, while the observed large-scale 
motions have not yet developed into galactic outflows. On 
the other hand, among 29 LBGs drawn from the bluest 
quartile of intrinsic UV colors, only two of them have their 
starburst ages estimated to be ~ 1 Gyr under the assumption 
of constant star formation (Shapley et al. 2001). However, 
there seems to be a relationship such that LBGs with 
intrinsically bluer rest-frame UV colors, after correction for 
IGM absorption, are older starburst systems (Shapley et al. 
2001). According to the results with our model disks with 
Mrf = 10* Mq at high redshift, we could suggest that the 
escape fraction can easily be ~0.5 in the LBGs, because the 
strong galactic outflows due to powerful starbursts can 
make chimneys allowing a significant number of ionizing 
photons to escape to the IGM. Many LBGs show evidence 
for large-scale galactic outflows (Pettini et al. 1998, 2001). 
However, Fernandez-Soto et al. (2003) do not find the trend 
of increasing escape fractions as a function of color of gal- 
axies. The low escape fractions observed in them are hard to 
explain with our models, unless they are all starburst sys- 
tems so young that the shells still trap the ionizing radiation. 

6.2. UV Background Radiation 

Our results suggest that a significant fraction of ionizing 
radiation escapes from high-redshift disks. The recent mea- 
surement by the Wilkinson Microwave Anisotropy Probe of a 
large electron scattering optical depth implies that a signifi- 
cant volume of the IGM was reionized by redshift z ~ 17 ± 5 
(Kogut et al. 2003; Spergel et al. 2003). Population III objects 
may account for such early reionization if efficient molecular 
hydrogen cooling occurs. However, the measured tempera- 
ture of the Lya forest at z ~ 3 rules out a single reionization 
model before z ~ 9 and suggests the reionization history is 
complex (Haiman & Holder 2003; Wyithe & Loeb 2003; Hui 
& Haiman 2003; Theuns et al. 2002). Therefore, we attempt 
to make an order-of-magnitude calculation of the contribu- 
tion of escaping UV photons from dwarf starburst galaxies at 
z > 5 to the UV background radiation field and thus to the 
reionization of the universe. 

We compute as a function of redshift the rest-frame 
comoving density of Lyman continuum photons escaping 
out of dwarf galactic halos. To do so, we first compute the 
formation rate of halos -Rform and the survival probability of 
the halos Psurv as a function of redshift, using the excursion 
set halo mass function for a hierarchical random Gaussian 
density field (Bond et al. 1991) and the theory of hierarchical 
merging probabilities (Lacey & Cole 1994). Refer to Lacey 
& Cole (1994) for the detailed derivation of TJform and Psurv- 
We use the spherical collapse model to define the critical 
density threshold for halo collapse and the CDM power 
spectrum given by Bardeen et al. (1986) with a scale- 
invariant initial condition to define the mass variance, as in 
Kitayama & Suto (1996). We employ a ACDM model with 
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fio = 0.37, = 0.63, Qb = 0.05, h = 0.7, and ag = 0.8, 
although /esc(0 and F^sc were computed with our standard 
model (SCDM+(/)i(;i) in the previous sections. Recall that 
the computation of -fesc is not sensitive to the change in 
cosmology nor dark-to-baryonic ratio within theoretical 
and observational uncertainties. 

We select two groups of halos. First we look at halos with 
virial temperature 71 > lO"* K [minimum halo mass 
^min(z) = A/4(z)], in which moderately efficient star forma- 
tion is possible because of hydrogen line cooling. Second, 
we consider halos with > 10^ K or > 50 km s~' 
[Mnin(z) = Af5(z)], in which >50% of the gas within the 
virial radius can cool in the reionized IGM (e.g., Thoul & 
Weinberg 1996). The maximum halo mass for the computa- 
tion is set to M,nax = 10" Mq,. We assume the gas cools and 
settles into disks promptly at the redshift of halo collapse 
because the cooling time of the halos is much smaller than 
the Hubble time at high redshift, ^cooi ^ ?h (Oh & Haiman 
2002). 

We then define the Lyman continuum photon back- 
ground density due to dwarf starburst galaxies with star 
formation timescale tsf as 



Puv 



clM 



dzfS{M, Zf, z)f^ 



^ Fesc(M,/^,z)6ph(M,/^) 

A J Tsf 



(18) 



where Mis the halo mass,/=i= is the star formation efficiency, 

assumed to be constant in M and z, Qph{M,ft.) is the total 
number of photons produced per unit gram of gas converted 
to stars, and S{M,Zf,z)dM dzf = S{M , tf, t)dM dtj- is the 
number density of halos that form with mass M ~ M + dM 
at time tf ~ tf + dtf and survive without destruction until a 
later time t (Lacey & Cole 1994). This can be expressed as 

S{M, tf, t)dMdtf = R[orm{M, t)dM dtf Psmv{M, tf, t) . 



(19) 

Equations (18) and (19) simplify to 

X i^esc(M,/*, z)eph(M,/*)J?form(M, z)dM (20) 

if we assume instantaneous starbursts, since most photons 
are produced within the first few Myr of the bursts. 
Equation (20) also reproduces equations (18) and (19) with 
multiple starbursts having tsf = 20 Myr because ^form and 
Psurv do not change over such a short timescale. 

At z > 5, we assume Fesc(M,/*,z) to be constant based 
on our results with Mj =10^ Mq-, at z = 3 and 8 with a star 
formation efficiency /* = 0.06. To span the range of possi- 
bilities, we consider total escape fractions i^esc = 0.2, 0.4, 
and 0.8, so that f^Fesc takes on values of 0.012, 0.024, and 
0.048. We also choose />i==0.1 and Fesc = 1.0, so 
/*-fesc = 0.1 as an upper limit. 

We show in Figure 14 the evolution of the UV back- 
ground radiation as a function of redshift for z > 5 under 
our assumptions, compared with the minimum photon pro- 
duction rate needed to fully ionize the universe (i.e., when 
the ionization fronts completely overlap). The minimum 
rate is defined as J/" = «hi(z = 0)/7rec, where mhi(z = 0) is 
the comoving mean hydrogen density of the expanding 
IGM and the volimie-averaged recombination time 
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Fig. 14. — UV background density (photons s^' comoving Mpc"') 
computed from two groups of halos (1) with !"„ > IC K (M > A/4; solid 
line) and (2) with T, > 10^ K (M > Ms; triple-dot-dashed line), based on 
the escape fractions F^^c in disks with Mj = 10* Mq at r = 3 and 8 studied 
above. The photon production rate is proportional to /*-Fcsc when we 
assume that both /* and fcsc are constant in z and Mj. It has to exceed the 
minimum photon production rate to fully ionize the imiverse if we assume 
clumping factors C = 30, 50, or 100 (dashed lines). 



ttec = [«Hi(z)<a;BC] with the ionized hydrogen clumping 
factor C (Madau et al. 1999; Gnedin & Ostriker 1997). This 
is based on the assumption that the clumping can be aver- 
aged over, and that the recombination time is very small, 
'rec^ ^h(z) at high redshift so that only photons emitted 
within time 7rec can actually ionize new material. We take 
values of C ranging from 30 to 100 on the basis of numerical 
simulations (Gnedin 2000) and analytical studies (Benson 
et al. 2001; Haiman, Abel, & Madau 2001). 

Figure 14 suggests that the Lyman continuum photon 
production rate from dwarf starburst galaxies are high 
enough to make a substantial contribution to reionize the 
IGM before z « 10 if /*f;sc ^0.024; that is,/* >0.06. We 
argued that /* >0.06 is plausible on the basis of gravita- 
tional instability of our model disks and also with the pres- 
ence of metals replenished from previous generations of 
stars. However, the formation of dwarf galaxies can be sup- 
pressed by mechanical feedback from other dwarf galaxies 
(Scannapieco, Thacker, & Davis 2001; Fujita et al. 2003). 
To study the process of reionization, we ultimately need 
three-dimensional cosmological simulations with radiative 
transfer, which can also resolve the formation of SN-driven 
bubbles and outflows. Our crude calculation still shows that 
dwarf starburst galaxies with i^esc^O.2 with a reasonable 
star formation efficiency,/* « 0.05-0.1, provide a substan- 
tial amount of UV photons to the IGM. 

7. CONCLUSION 

We used numerical models to study the influence of shells 
and galactic outflows on the escape of ionizing radiation 
from dwarf galactic disks with starburst clusters. From our 
models, we can draw the following conclusions: 

1. The shells of ISM swept up by isothermal shocks 
around young superbubbles can trap ionizing radiation very 
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effectively until the bubbles start to accelerate, causing the 
shells to fragment. The acceleration of the bubbles begins 
when the bubbles expand to a few scale heights above the 
disks. 

2. The low escape fractions of ionizing photons observed 
in local dwarf galaxies less than 0.06 can be explained if they 
are observed before blowout. However, the observed fesdt) 
may not reflect the total fraction of photons escaping from 
the disks over the lifetime of the starbursts F^qsc • 

3. Galactic outflows blow out to form holes that allow 
ionizing radiation to escape directly to the IGM. The escape 
fractions are ^esc ~ 0.2 in the high-redshift disks, in which 
only a negligible amount of photons can escape without 
bubble dynamics. This is because the gas distributions are 
highly stratified in the strong gravitational potentials of 
dense, high-redshift halos, allowing quick acceleration of 
the bubbles. 

4. With =0.05-0.1 and Fesc>0-2, dwarf starburst 
galaxies may have contributed significantly to the UV back- 
ground radiation and so to the reionization of the IGM 
before z 10. 

The assumptions used in modeling dwarf galaxy disks 
and star formation distributions and histories can influence 
the predicted escape of ionizing radiation in several ways: 

1. The vertical scale height H influences the escape of 
ionizing photons because the recombination measure is 
Sz oc Njj^. Our thick undisturbed disks suppress escape 
more than our thin undisturbed disks. Galactic outflows 
enhance this effect because bubbles blow out of thin disks 
earlier than thick disks: ?bi oc Z/"*/^. The time of blowout 
determines the transition from the trapping of photons by 
shells to their escape through chimneys in galactic disks. 

2. In undisturbed disks, the peak photon luminosity pri- 
marily determines the escape fractions. The instantaneous 
escape fraction /esc(f) follows the behavior of A^ph(?), pre- 
dicted by a star formation history (e.g., instantaneous burst 
or multiple bursts). In disks with supershells and galactic 
outflows, fesc{t) strongly depends on the evolution of 
bubbles rather than on the star formation history. The 
escape fraction with a different star formation history from 
the ones we used in this study can be inferred on the basis of 
our results by computing and comparing its peak photon 
luminosity and expected mechanical luminosity. 

3. The star formation efficiency /* we used in this study is 
for a single starburst cluster in a disk, not for the entire disk. 
When a fixed total star formation efficiency is given in the 
disk, the distribution of ionizing sources strongly influences 
the escape of ionizing photons because the photon luminos- 
ity from each ionizing source must exceed a minimum pho- 
ton luminosity Np^^^i^ for its Stromgren sphere to break out 
of the disk. A single ionizing source gives the most favorable 
condition for the photons to escape, compared with numer- 
ous OB associations scattered around the disk in space and 
in time, as considered in Dove & Shull (1994), Dove et al. 
(2000), Wood & Loeb (2000), and Ciardi et al. (2002). 

Our models have several limitations: 
1. Shell formation is not resolved even with resolutions 
as low as a few tenths of a parsec. Therefore, when we com- 



pute the escape of ionizing radiation in the simulations, we 
use the theoretical peak density to compute the recombina- 
tion measure H along rays that encounter dense shells rather 
than holes produced by the Rayleigh-Taylor instability. 
This correction was necessary only in low-redshift disks 
with low hydrogen density. 

2. The fragmentation of the shells is resolved qualita- 
tively: the hot gas accelerates beyond the cold swept-up 
shells. The opening angles of chimneys cleared out by the 
Rayleigh-Taylor instability remain 30°-50° at the time of 
blowouts regardless of resolution. However, the quantita- 
tive details of the fragmentation, including the number and 
positions of individual fragments, are determined by our 
assumption of cylindrical symmetry, by numerical 
resolution, and by the lack of magnetic fields in our 
calculations. 

3. We modeled only a single starburst cluster placed at 

the center of a rotationally supported disk. However, we 
argue that a starburst cluster placed at the center yields a 
lower limit for the escape fraction, which remains as a lower 
limit in a more realistic galaxy in which multiple star clusters 
of the same or smaller photon limiinosities are scattered 
around the disk. 

4. Our models include only a single-phase, smoothly 
stratified ISM. We argue that this yields a lower limit for the 
actual escape fraction for ionizing radiation. 

5. We assume weU-formed, rotationally supported disks, 
even at high redshift. A more turbulent ISM may form 
because of stellar feedback and cooling, perhaps even before 
the formation of a disk. Bubble blowout through such a 
medium is likely to be easier than through a well-formed 
disk. However, to study the escape fractions at high redshift 
and the reionization of the IGM, it will ultimately be 
necessary to follow the collapse of gas and solve the radia- 
tion transfer problem in three-dimensional cosmological 
simulations. 
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